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Dispositif d'analyse photopyroelectrlque. 

fej) Dispositif d'anaiyse photopyroelectrique, caracterise 
en ce qu'il comports des moyens (7) de calcul dans les- 
quels est memorises une model isation bidimensionnelle 
axisyrnetrlque des transferts thermiques qui Intervlennent 
dans r ensemble multicouche constitue par Techantillon (2) 
et le capteur pyroelectrique (4). 
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La pr^sente invention est relative a un dispositif 
pour l f analyse thermique photopyroeiectrique de mat^riaux 
structures en couches. 

Les m6thodes photothermiques , de d6veloppement 
relativement r6cent, sont particuli&rement bien adapt£es 
pour 1' analyse de la structure et des propri6t6s 
thermophysiques de materiaux minces en couches. Elles 
trouvent en particulier applications dans le cas des 
mat£riaux solides, mais peuvent aussi etre employees pour 
la caract6risation de liquides ou de gaz. 

Elles consistent h soumettre un 6chantillon St 
analyser ^ un rayonnement excitateur module, par exemple a 
l'aide d'un hacheur m6canique, et suivre Involution 
d'un signal photothermique dependant de la temperature 
d'une zone de 1 9 echantillon. L'amplitude du signal 
photothermique mesure et son dephasage par rapport au 
rayonnement excitateur permettent d' avoir acc6s Si certains 
parametres de 1 ' echantillon, tels que son absorptivity, sa 
diffusivite, sa structure et ses d6fauts internes f etc. 

Parmi les techniques utilis6es jusqu'^ present 
pour mesurer un signal photothermique, on peut citer : 
. la radiometrie infrarouge ou l'on capte les variations 
d'emittance de la zone chauffee i l'aide d'une chaine de 
detection radiometrique infrarouge (capteur ponctuel ou 
camera de thermographie ) ; des applications particulidres 
de la radiometrie photothermique ont notamment ete 
decrites dans les brevets frangais deiivres au nom de la 
Demanderesse et publics sous les numeros 2 593 917 
et 2 663 745 ; 

. l'effet mirage, ou on utilise un faisceau laser "sonde", 
traversant la couche d'air au voisinage de la surface, 
qui se trouve devie par le gradient d'indice de 
refraction produit par l'616vation de temperature. 
Beaucoup plus sensible que la mfrthode radiometrique, ce 
proc6d6 de detection presente par ailleurs quelques 
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inconvenients, -notamment celui du reglage de la distance 
du faisceau sonde par rapport a 1 ' echantillon ; 
les capteurs utilisant des materiaux pyro<§lectriques en 
films minces, egalement de grande sensibilite et en.meme 
temps de faible cout. lis sont d'un emploi aise dans le 
cas de materiaux non-solides (liquides, pates, gels 

La detection pyroelectrique de plusieurs ordres de 
grandeurs plus sensible que la radiometric infrarouge, 
permet d'obtenir des signaux exploitables avec des 
puissances d 'excitation environ 100 fois plus petites . 
L' analyse photopyroelectrique est done particulierement 
avantageuse pour 1' analyse de materiaux thermiquement 
instables, tels que des liquides qui peuvent etre l'objet 
de phenomenes de convection. 

En outre, le cout des moyens d' excitation, ainsi 
que leur taille, se trouvent considerablement r6duits. 

Ces avantages ouvrent des perspectives nouvelles 
d' applications pour 1' analyse photothermique . Les 
techniques pyroelectriques permettent d'envisager des 
applications de 1' analyse photothermique au domaine 
medical (sedimentation, coagulation du sang), a la 
biologie (developpement de couches cellulaires) , a la 
mineralogie (etudes de granulometr ie ) , a la chimie (etudes 
de changements de phase). 

Jusqu'a present, cependant, les techniques 
photopyroelectriques ne permettaient de realiser que des 
analyses empiriques, par comparaison avec des donnees 
experimentales preetablies. Les mesures realisees 6taient 
confrontees a des abaques de reseaux de courbes fonction 
des valeurs du ou des parametres physiques a determiner. 

L' invention propose quant a elle un dispositif 
permettant des analyses quantitatives d 'echantillons 
structures en couches minces, sans qu'il soit necessaire 
d'etablir, avant toute mesure et pour chaque type 



2715226 



3 

d'£chantillon, des r^seaux de courbes exp^rimentales de 
comparaison. A cet effet, le dispositif selon 1' invention 
met en oeuvre une mode lis at ion adaptee des comportements 
thermiques des 6chantillons . 
5 II a d6j§i ete propose dans les techniques 

photopyroeiectriques de mettre en oeuvre des modeiisations 
unidimensionnelles des transferts thermiques se produisant 
dans les echantillons analyses. On pourra ^ cet egard 
avantageusement se referer aux publications suivantes : 
10 * A. MANDELIS, M. M. ZVER, "THEORY OF PHOTOPYROELECTRIC 
SPECTROSCOPY OF SOLIDS" - J. Appl . Phys. 51, 4421 
(1985). 

* M. CHIRTOC, G. MIHAILESCU, Phys. Rev. B40, 9606, (1989), 
Ces mod61isations unidimensionnelles sont 
15 neanmoins insuf f isantes pour rendre compte des trans ferts 
thermiques intervenant, lorsque l'on utilise un faisceau 
d' excitation cylindrique de rayon faible. 

L' invention propose quant ^ elle de mettre en 
oeuvre une modeiisation bidimensionnelle axisymetrique des 
20 trans ferts thermiques intervenant dans les echantillons 
multicouches analyses. 

Une telle modeiisation est d6ji connue et a dej& 
ete utilisee notamment en radiom6trie photothermique. On 
pourra h cet egard avantageusement se r6f6rer a l'ouvrage 
25 de th^se : "Radiometric photothermique appliquSe a la 
caracterisation et au controle du traitement thermique des 
metaux" F- POTIER, 30 novembre 1989 - UNI VERS I TE DE REIMS 
CHAMPAGNE - ARDENNE - Sp6cialit6 energ6tique. 

A ce jour, une telle modeiisation n'avait pas ete 
30 utilisee pour 1' analyse thermique photopyroeiectrique. 

Elle s'avdre n6anmoins trfes bien corr616e avec 
1 ' experimentation . 

En outre, une telle modeiisation permet egalement 
de prendre en compte les trans ferts thermiques qui 
35 interviennent entre le capteur pyroeiectrique et 
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1' echantillon analyst, aihsi que dans le capteur 
pyroelectrique lui-meme, 1 'echantillon et le capteur 
pyroelectrique etant modelises comme un ensemble 
multicouche . 

L' invention propose un dispositif d' analyse 
photothermique d'un materiau structure en couches, pour la 
determination d'au moins un parametre physique d'une de 
ces couches, tel que son epaisseur, sa diffusivite, son 
absorption optique, comportant: 

- une source pour 1 'emission d'un signal lumineux 
d' excitation en direction d'un echantillon du mat6riau a 
analyser, la direction d' incidence de ce signal par 
rapport a 1 'echantillon etant perpendiculaire aux plans 
des couches de celui-ci ; 

- des moyens pour la modulation en amplitude de ce signal 

d' excitation, 

- un capteur pyroelectrique en contact thermique avec une 
face de 1 'echantillon, cette face etant parallele aux 
plans desdites couches, 

- des moyens pour mesurer un signal photothermique de 
tension aux bornes de ce capteur pyroelectrique, 

caracterise en ce qu'il comporte en outre : 

- .des moyens pour determiner au moins une valeur numerique 
representative de 1' amplitude dudit signal 
photothermique et/ou au moins une valeur numerique 
representative du dephasage entre ce signal 
photothermique et le signal d' excitation ; 

- des moyens de calcul dans lesquels est memorisee une 
modelisation bidimensionnelle axisymetrique des 
trans ferts thermiques qui interviennent dans 1' ensemble 
multicouche constitue par 1 'echantillon et le capteur, 
ces moyens de calcul permettant de determiner la valeur 
du parametre physique correspondant au signal 
photothermique de tension mesure. 

D'autres caracteristiques et avantages.de 
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1' invention ressortiront encore de la description qui 
suit* Cette description est purement illustrative et non 
limitative. Elle doit etre lue en regard des dessins 
annexes sur lesquels : 
5 - la figure 1 represente schematiquement un 

dispositif pour la mise en oeuvre de 1' invention ; 

- les figures 2 et 3 representent schematiquement 
deux configurations d' analyse possibles ; 

- la figure 4 illustre la mod61isation thermique 
10 mise en oeuvre par le dispositif conforme £ 1' invention ; 

- la figure 5 illustre une configuration d' analyse 
mise en oeuvre pour verifier la correlation entre la 
modeiisation propos^e par 1 ' invention et les r^sultats 
experiment aux ; 

15 - les figures 6 et 7 sont des graphes sur lesquels 

ont ete portees, en fonction du paramfetre X represents sur 
la figure 5, d'une part, 1' amplitude et la phase du signal 
photothermique de tension releve avec le montage de la 
figure 5 r et, d' autre part, les courbes theoriques 

20 d' amplitude et de phase theoriques du signal de tension 
calcuie h partir de la modeiisation mise en oeuvre par le 
dispositif de 1' invention. 

Le dispositif represents sur la figure 1 comporte 
une source 1 d' excitation d'un milliwatt de puissance, 

25 emettant un faisceau F de lumi&re monochromatique sur un 
echantillon 2 Si analyser. En sortie de cette source 1, le 
faisceau F est module par des moyens qui peuvent etre 
constitues par un hacheur mecanique class ique, ou encore 
ainsi qu'on l'a repr6sente sur la figure 1 par un 

30 modulateur eiectro-optique command^ par une unite 3 lui 
deiivrant un signal de modulation. D'autres sources 
lumineuses d' excitation sont bien entendu envisageables . 
En particulier, une excitation par une simple LED suffit. 
D'autres types de modulations sont 6galement possibles, 

35 telles qu'en particulier une modulation par techniques 
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acousto-optiques . 

Sur une face de 1 ' 6chantillon 2, est monte un 



polyf luorure de vinylid6ne (PVDF). Le signal de tension 
aux bornes de ce capteur 4 est envoy£ sur un 



a un amplif icateur 6 i detection synchrone du type 
Stanford SR-850. 

Les tensions 6chantillonn6es par 1 'amplif icateur 6 
sont envoy£es sur les moyens de calcul 7 h 1' entree 
desquels ils sont num£ris£s. Les moyens de calcul 7 
"traitent 1' amplitude et la phase des signaux regus pour 
determiner le param&tre physique de 1 '^chantillon que l'on 
cherche h mesurer. 

On se r6f6re maintenant plus particuli&rement aux 
figures 2 et 3 sur lesquelles ont 6t6 illustrees deux 
configurations d' analyse possibles, 

Dans la configuration de la figure 2, 
l'6chantillon multicouche 2 re<joit directement le faisceau 
F module, le capteur pyro61ectrique 4 etant dispose sur la 
face de 1 '6chantillon 2 oppos6e a celle qui re^oit le flux 
lumineux. 

Dans une autre configuration possible, appeiee 
"configuration inverse", c'est le capteur pyroeiectrique 4 
qui regoit le faisceau F. Sa face dirig6e en regard de la 
source 1 porte une couche noire 8 absorbant l'6nergie 
calorifique du faisceau F, 

Pour une description plus d6taill£e des capteurs 
pyroeiectriques et des montages qui leur sont associ6s, on 
pourra avantageusement se r£f6rer aux articles : 
- "A versatile inverse photopyroelectric (IPPE) technique 
and instrument for real time observation of the 
condensation of water vapor in the atmosphere" - Mihai 



capteur pyroeiectrique 4. Ce capteur pyroeiectrique 4 est 
par exemple du type classique & feuille mince de 
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CHIRTOC, Dane BICANIC, Valer TOSA - Rev. Sci. Inst rum. 
62 (9), September 1991 pp. 2257-2261 ; 
- "Pyroelectric Sensors for the Photothermal Analysis of 
condensed Phases - H. COUFAL and A. MANDELIS, 
5 Ferroelectrics 1991, vol. 118, pp. 379-409. 

Classiquement, le capteur 4 est compost d'un film 
pyroeiectrique plan de surface A et de largeur L, 
m4tallisd sur ses deux faces (depots de faible epaisseur). 
Ce film est caracterise par une polarisation eiectrique 
10 spontanee P, fonction de la temperature T. II resulte de 
cette polarisation une apparition en surface de charges 
f ictives Q dont la density surf acique a est egale §L la 
compos ante normale de la polarisation Pn- Ces charges, 

liees, creent un champ eiectrique interne qui, par le 
15 deplacement des charges libres qu'il entraine, aboutit h 
la neutralisation de la charge superf icielle . Les 
perturbations d'ordre thermique modifient cet 6quilibre de 
surface en agissant sur la polarisation. Ces variations 
thermiques de la polarisation sont caract6ris6es , h. une 
20 temperature donnee, par le coefficient pyroeiectrique p 
(en Cm-2K-i). Toute variation de charges li6es dans le 
dieiectrique, consecutive 2l une variation de la 
[temperature, peut ainsi etre d6tect6e en mesurant la 
difference de potentiel qui s'etablit entre les armatures 
25 metalliques . 

La tension aux bornes du capteur pyro61ectrique 
[jest donnee par : 

[it J 

30 avec 

■6(0 «fc |T(r-x)dx 
ou e r est la permittivite relative caracteristique du 
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dieiectrique (6gale environ k 12 pour le PVDF) . 

On se ref&re maintenant k la figure 4. 
Confona^ment k l r invention, on considfere que 1 ' echantillon 
2 et le capteur 4 en contact avec celui-ci constituent un 
5 element a symetrie axiale a trans ferts thermiques 
bidimensionnels. L'axe de symetrie est l'axe du faisceau 
d' excitation F, perpendiculaire aux plans des couches. 

La modeiisation bidimensionnelle et axisymetrique 
des transferts thermiques qui s'effectuent au sein de 

10 1' ensemble du syst&me echantillon/capteur permet de 
determiner le champ de temperature dans cette structure") 
multicouches et d'6tablir une relation de calcul entre le 
signal de tension photopyroeiectrique et le signal 
d' excitation envoye sur 1 'echantillon. _J 

15 l' echantillon que l'on veut etudier est suppose en 

bon contact thermique sur une surface S de 1' armature du 
detecteur. L ' absorption lumineuse dans ce systSme 
echantillon/detecteur et la conversion supposee non 
radiative de l'energie lumineuse en chaleur , augmente la 

20 temperature du film pyroeiectrique (dans une region que 
l'on suppose limit6e k la surface S). Si 1' excitation est 
moduiee sinusoldalement, la variation de temperature qui 
en resulte, entraine k son tour une variation periodique 
de la charge accumuiee dans le pyroeiectrique , dephas6e 

25 d'un angle $ P ar rapport au flux excitateur. 

La modeiisation bi-dimensionnelle et axisymetrique 
des transferts thermiques qui s'effectuent au sein de 
1' ensemble du systeme 6chantillon/d6tecteur , a permis dans 
un premier temps, de determiner le champ de temperature 

30 dans cette structure multicouches et de calculer ensuite 
le signal photopyroeiectrique. 

L' expression du flux excitateur est : 

I|(M) ^I (r) l±|Jifi 

35 Lorsque 1' excitation est un laser £ profil 
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gaussien, 1 'eclairement est de la forme : 

i(r>- r 9 t 43 

ou 10 est la density de puissance au centre du faisceau 
5 laser, ro le rayon d' excitation pour lequel la valeur de 
<(>q est divisee par e (10 = P / **q2). Le faisceau est 
normal a la surface ; il faut tenir compte des flux 
radiatifs reflechis Ii~ et transmis I ± + se propageant 

parallelement a l'axe du cylindre, ceci dans chacune des 
10 couches du systeme. Il faut alors resoudre le systeme de 
deux equations issues de 1' equation de transfert dans les 
deux directions particulieres . Si i est l'indice de la 
couche et pi le coefficient d' absorption optique donne 
relativement a cette couche, on a alors : 



15 



20 



25 



dx 

j r ' (r *l . Pir,(r.r) -0 



dont les solutions sont : 



If(r.z)-Pie*(^L-P*« 
Ii(r,z)«Mje-(riLp* 



On n6glige ici la diffusion de la lumiere. Les 
30 conditions aux limites en z nous permettent de calculer 
les coefficients Pi et Mj_, et par suite, d'en d6duire 

1' expression des sources de chaleur g^neree au sein de 
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chacune des couches i : 

C,(i\z)«- dhr lf(rjt) 
E|(f\i)« - div Ij(r.z) 

Les flux radiatifs n'ayant pas ici de composantes 
radiales, ces expressions deviennent dans le cas d'une 
source moduiee de pulsation a) : 



10 CKr t z)-l^,P,e-[3r^J«t 



On suppose que l f absorption de l'energie est 
15 totale et que celle-ci est integralement convertie en 
chaleur. On suppose realiser les conditions d'adiabaticite 
des parois laterales. La resolution de 1' equation de la 
chaleur dans chacune des couches nous permet alors de 
determiner compietement le champ de temperature dans 
20 1 'echantillon et dans le pyroeiectrique. Cette equation 
relativement ^ une couche i a pour expression : 



25 ou 

- Ki est la conductivite thermique de la couche i, 

- Ti la temperature de la couche i, 

- ri la masse volumique, 

- pi la chaleur specifique, 

30 - Gi et Ei les termes de sources de chaleur g6n6r6e dans 
la couche i. 

La temperature s'exprime alors comme la somme de 
trois composantes : 

T|(r^4) • T„b+T c (r*HT.t r,z)ei" 
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ou 



10 



Tamb est la temperature ambiante, 
T c (r,z) la composante continue de la temperature, 
T a (r,Z) la composante alternative en tout point de la 
couche i. 

Pour la determination du signal 
photopyroeiectrique, seule la composante dynamique de la 
temperature est & prendre en compte. Elle a pour 
expression : 

Trfr,z>.£ J^r)[A|..e«»-^Bi..e-«'**.Eue»' s -G Lt e ^| 



••0 

15 ou 

r c 6tant le rayon de 1 'echantillon, 
bn est le n-idme zero de Ji, 

ai est la dif f u'sivite' du milieu dans la couche i, 

P|C| 

25 et enfin, Pi est le coefficient d' absorption optique de la 
couche i. 

Les constantes Ai, n sont determinees par les 

conditions aux limites sur z, & savoir s 

- conditions mixtes et coefficients d'echange en 
30 face avant et arriere, 

- conservation du flux et 6galit6 des temperatures 
aux interfaces. 



BNSDOC1D: <FR 



271S226A1J_> 



2715226 



12 

Les termes Ei /n G±,n sont fonctions entre autre 
de coefficients I n obtenus en d6composant les termes de 
sources sur la base des fonctions de Bessel telles 
que Jo : 

p. 



Pi 



Mil. 
Pil. 



10 avec 



I. -II- 

JS(X.r t ) 



15 



20 



Si on choisit pour i l'indice correspondant £ la 
couche du pyro^lectrique, on peut alors en d6duire 
1 'expression de la tension pyro61ectrique s 



VppeM 



XjQ _ Yjp 
<*u> Pi 



hQ^c) \X in Y 



+ 2 £ , 

n=l Kfc l°in Pi 



A(Cu) d** 



25 avec 



30 



X in = Ain(e-^ - e-°^') + B in (e - e 



Yin - E in (e - P* - e - P*>) + G in (e P*-. - c P*) 
q, le rapport des surfaces echantillon/d6tecteur, 
A(to), la fonction de trans fert du r6seau 61ectrique fonn6 
du d^tecteur connect^ a un pr6amplif icateur convertisseur 
courant-tens ion . 

II s'agit ici de quant ites complexes. On peut 
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alors ecrire cette expression en coordonn6es polaires et 
en tirer les valeurs de 1' amplitude et de la phase du 
signal photopyroeiectrique. 

Les moyens de calcul 7 memorisent cette 
model isat ion et sont programmes de fagon h resoudre cette 
equation pour determiner, £ partir de 1' amplitude et de la 
phase qu'ils regoivent, le paramStre physique que l'on 
cherche a mesurer. 

Dans une variante possible, les moyens de calcul 
determinent des courbes theoriques donnant la phase ou 
1' amplitude en fonction du paramStre que l'on cherche i 
mesurer, et l'on determine ce paramdtre & partir de la 
phase ou de 1' amplitude mesurees, en utilisant ces courbes 
comme abaque* 

On notera que la modeiisation propos6e permet de 
determiner la tension th6orique aux bornes du capteur 
pyroeiectrique & partir de la connaissance, d'une part, de 
la puissance et de la modulation du faisceau d' excitation 
et, d' autre part, de la diffusivite thermique, de 
1' absorption et de 1'epaisseur de chaque couche de 
1' ensemble echantillon/capteur. 

Le dispositif conforme & 1' invention est 
avantageusement utilise pour le suivi des changements 
d'une structure physique et notamment pour le suivi de 
phenomenes de sedimentation ou de croissance de couche* 
Dans ce cas, pr6f 6rentiellement mais non n6cessairement 
l'echantillon analyse est modeiise comme un stratifie 
bicouche. La longueur d'onde du faisceau d' excitation est 
choisie de fagon a differencier les deux couches 
adjacentes. L' invention est en particulier avantageusement 
mise en oeuvre pour le suivi de la sedimentation ou de la 
coagulation d ' echantillons de sang* 

Dans le cas du suivi d'une sedimentation 
erythrocytaire, on peut en premiere approximation 
consid6rer 1 '6chantillon sanguin comme un stratifie bi- 
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couche, constitu£ d'une couche de sang pratiquement opaque 
h 1 'excitation et d'une couche de plasma semi-transparente 
a la longueur d'onde du faisceau d' excitation et dont 
l'£paisseur croit ayec le temps. Connaissant notamment la 
conductivity thermique du sang et du plasma (resp. 0,58 et 
0,56 W/m/K), leur masse volumique (resp. 1026 et 1108 
kg/m3), leur chaleur sp6cifique (4180 J/kg/K k la fois 
pour le plasma et le sang) , il est possible de determiner 
k tout instant, en fonction des mesures de 1' amplitude et 
de la phase du signal photothermique, l'epaisseur de la 
couche de plasma. L'epaisseur totale de 1 'ensemble 
constitue par 1 .' echantillon k analyser et le capteur 
pyroeiectrique est de l'ordre de 1 mm. 

Pour 1' analyse d'un suivi de coagulation , on 
travaille pr6f erentiellement en configuration inverse de 
fa£on que 1 '6chantillon ne soit pas perturbe optiquement . 

L' invention trouve egalement avantageusement 
d'autres applications. En particulier, elle peut etre 
utilisee pour la detection d' apparitions de bacteries ou 
pour le suivi de la croissance de cellules, phenom&nes 
dont les techniques photothermiques classiques ne 
permettent pas de rendre compte. 

On se refSre maintenant & la figure 5. On a 
represent^ sur cette figure la configuration d' analyse 
avec laquelle il a 6t6 verifie que la mod61isation mise en 
oeuvre par le dispositif conforme k 1' invention etait 
convenablement correiee avec 1' experimentation. 

L' echantillon 2 est un echantillon d'eau dans 
lequel est dispos6e une feuille 9 d ' aluminium k surface 
noircie qu'il est possible d'6carter ou de rapprocher du 
capteur 4 . On a r6f 6renc6 par X la distance entre le 
capteur 4 et la surface 8. 

L ' amplitude V et la phase 0 du signal de tension 
mesure sur le capteur 4 ont 6te portees pour plusieurs 
valeurs de la distance X sur les graphes des figures 6 et 
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7 . On a egalement port6 sur ces graphes les courbes Vc et 
0c d' amplitude et de phase calcul6es k l'aide de la 
mod61isation mise en oeuvre par le dispositif conforme & 
1' invention. Ces graphes montrent que la mod£lisation 
5 propos^e permet ef fectivement d' avoir accfes par le calcul 
5t des valeurs bien corr616es aux valeurs mesur6es 
experiment alement . 

On a d6crit 1' invention dans le cas oh le faisceau 
lumineux envoy6 sur 1 'Schantillon est un faisceau moduli 

10 sinusoidalement • Bien entendu, 1' invention s' applique 
6galement dans le cas ou 1 'excitation de 1 '6chantillon est 
du type pseudo-al6atoire* On se r£f6rera 2k cet €gard 
avantageusement au brevet frangais de la Demanderesse 
public sous le num6ro 2 663 745. 

15 On notera qu'un des avantages de 1' invention est 

de permettre une analyse £ une frequence de modulation 
inf6rieure ct 1 Hz, ce qui permet de travailler avec des 
£paisseurs d'^chantillons allant jusqu'au mm. 
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REVENPICATIOHS 

1. Dispositif d' analyse photothermique d'un 
materiau structure en couches, pour la determination d'au 
moins un paramStre physique d'une de ces couches, tel que 
5 son 6paisseur, sa diffusivite, son absorption optique, 
comport ant : 

- une source (1) pour 1 'Emission d'un signal lumineux 
d 'excitation en direction d'un 6chantillon (2) du 
materiau a analyser, la direction d' incidence de ce 

10 signal par rapport 2l l'6chantillon etant perpendiculaire 
aux plans des couches de celui-ci ; 

- des moyens (3) pour la modulation en amplitude de ce 
signal d' excitation, 

- un capteur (4) pyroeiectrique en contact thermique avec 
15 une face de l'6chantillon, cette face etant par allele 

aux plans desdites couches, 

- des moyens (5, 6) pour mesurer un signal photothermique 
de tension aux bornes de ce capteur pyroeiectrique, 

caracterise en ce qu'il comporte en outre: 
20 - des moyens pour determiner au moins une valeur numfirique 
representative de 1' amplitude dudit signal 
photothermique et/ou au moins une valeur numerique 
representative du dephasage entre ce signal 
photothermique et le signal d' excitation ; 
25 - des moyens (7) de calcul dans lesquels est memorise© une 
modeiisation bidimensionnelle axisymetrique des 
transferts thermiques qui interviennent dans 1' ensemble 
multicouche const itu6 par 1 '6chantillon (2) et le 
capteur (4), ces moyens de calcul permettant de 
30 determiner la valeur du param&tre physique correspondent 
au signal photothermique de tension mesur6- 

2, Dispositif selon la revendication 1, 
caracterise en ce que les moyens (5, 6) pour mesurer le 
signal photothermique de tension et pour determiner une 
35 valeur num6rique representative de son amplitude et/ou de 
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son d£phasage comportent des moyens pour m6moriser les 
valeurs num6riques correspondant k des acquisitions de 
signal photothermiques r£alis6es & intervalles de temps 
r6guliers, de sorte qu f il est possible de r6aliser un 
suivi d'£chantillons evolutifs. 

3. Dispositif selon l'une des revendications 
pr6c6dentes, caract£ris6 en ce que les moyens de calcul 
mod61isent selon un flux gaussien le flux du signal 
d' excitation. 

4. Dispositif selon l'une des revendications 
pr£c6dentes, caract6ris6 en ce que la modeiisation 
m6moris6e dans les moyens de calcul prend en compte les 
flux radiatifs transmis et r£fl6chis dans chacune des 
couches de 1' ensemble const itu£ par 1 '6chantillon et le 
capteur. 

5. Dispositif selon l'une des revendications 
pr6c£dentes, pour l r analyse d'une sedimentation ou d'une 
croissance de couche. 

6. Dispositif selon l'une des revendications 
pr£c6dentes, caractfiris6 en ce que les moyens de calcul 
modSlisent l'6chantillon comme un stratififi bicouche. 

7. Dispositif selon l'une des revendications 5 ou 
6, pour 1' analyse de la sedimentation d'un echantillon de 
sang. 

8. Dispositif selon l'une des revendications 5 ou 
6, pour 1' analyse de la coagulation d'un echantillon de 
sang. 

9. Dispositif selon l'une des revendications 
pr£c6dentes, caracterise en ce que le signal d' excitation 
est moduli & une frequence de l'ordre de ou inferieure £ 
1 Hz. 
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The present invention relates to. an arrangement for the 
thermal photopyroelectric analysis of materials with a 
layered structure. 

5 Photothermal methods, which are a fairly recent 

development, are particularly suitable for the analysis of 
the structure and thermophysical properties of thin layer 
materials. They can be used in particular in the case of 
solid materials, but can also be employed for the • 
10 characterisation of liquids or gases. 

The methods consist in exposing a sample to be analysed to 
a modulated excitation radiation, for example by means of a 
mechanical chopper, and following the evolution of a 

15 photothermal signal depending on the temperature of a 
region of the sample. The amplitude of the measured 
photothermal signal and its phase difference with respect 
to the excitation radiation enable certain parameters of 
the sample to be evaluated, such as it absorptivity, its 

20 thermal dif fusivity, its structure and its internal 
defects, etc. 

Among the techniques used up to now for measuring a 
photothermal signal, the following may be mentioned: 

2 5 - Infrared radiometry, in which the variations in 

emissivity of the' heated zone are detected by means of 
an array of infrared radiometric detection devices 
(pinpoint sensor or thermographic camera) ; specific 
applications of photothermal radiometry have, in 

3 0 particular been described in the French patents granted 

in the name of the applicant and published under nos . 
2 593 917 and 2 663 745; 

The mirage effect, in which a "probe" laser beam, is 
used that traverses the layer of air in the vicinity of 



2 

the surface, and which is deflected by the refractive 
index gradient produced by the rise in temperature. 
This detection method, which is much more sensitive 
than the radiometric method, has certain disadvantages, 
5 however, particularly as regards the regulation of the 

distance of the probe beam with respect to the sample; 
Sensors employing thin film pyroelectric materials, 
which are likewise extremely sensitive and are at the 
same time inexpensive. These may conveniently be used 
10 in the case of non-solid materials (liquids, pastes, 

gels, etc . ) . 

Pyroelectric detection of several orders of magnitude more 
sensitive than infrared radiometry enables useful signals 
15 to be obtained using excitation strengths that are 
approximately 100 times smaller. Photopyroelectric 
analysis is therefore particularly advantageous for 
analysing thermally unstable materials such as liquids that 
may exhibit convection phenomena. 

20 

Furthermore, the cost of the excitation means as well as 
their size are considerably reduced. 

These advantages open up new areas of applications for 
2 5 photothermal analysis. Pyroelectric techniques enable 

applications of photothermal analysis to be considered in 
the medical field (sedimentation, blood coagulation) , in 
biology (development of cellular layers) , in mineralogy 
c (granulometry studies) and in chemistry (phase change 
30 studies) . 

Up to now, however, photopyroelectric techniques have 
allowed only empirical analyses to be performed, by- 
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comparison with predetermined experimental data. The 
measurements that were made were compared with nomograms of 
families of curves representing values of the physical 
parameter or parameters to be determined. 

5 

The invention per se provides an arrangement enabling 
quantitative analyses of samples with a thin layer 
structure to be carried out without having to establish 
families of experimental comparison curves before every 
10 measurement and for each type of sample. To this end the 

arrangement according to the invention employs a model that 
is adapted to the thermal behaviour of the samples.. 

It has already been proposed in photopyroelectric 
15 techniques to employ one -dimensional models of the thermal 
transfers occurring in the analysed samples. In this 
connection reference may advantageously be made to the 
following publications: 

* A. MANDELIS, M.M. ZVER, "THEORY OF PHOTOPYROELECTRIC ' 

20 SPECTROSCOPY OF SOLIDS" - J. Appl . Phys . 57, 4421 (1985). 

* M. CHIRTOC, G. MIHAILESCU, Phys. RevB40, 9606 (1989). 

These one-dimensional models are, however, insufficient to 
take account of thermal transfers that occur when a 
25 cylindrical excitation beam of small radius is used. 

The invention per se proposes to use an axially symmetrical 
two-dimensional model of the thermal transfers taking place 
in the analysed multilayer samples. 

30 

Such a model is already known and has already been used in 
particular in photothermal radiometry. Reference may 
advantageously be made in this connection to the thesis: 



4 

"Photothermal Radiometry applied to the. characterisation 
and monitoring of the thermal treatment of metals" 
F. POTIER, 30 November 1989 - UNIVERSITY OF RHEIMS, 
CHAMPAGNE -ARDENNES, Energetics Department. 

5 

Up to now such a model has not been used for thermal 
photopyroelectric analysis. 

Nevertheless, it is found that there is an extremely good 
10 correlation between this model and experimental values. 

Furthermore, such a model also enables the thermal 
transfers that take place between the pyroelectric sensor 
and the analysed sample as well as in the pyroelectric 
15 sensor itself to be taken into account, the sample and the 
pyroelectric sensor being modelled as a multilayer 
assembly. 

The invention proposes an arrangement for the photothermal 
analysis of a material with a layered structure, in order 
to determine at least one physical parameter of one of 
these layers, such as its thickness, its thermal 
diffusivity or its optical absorption, such arrangement 
comprising : 

a source for the emission of a luminous excitation 
signal in the direction of a sample of the material to 
be analysed, the direction of incidence of this signal 
with respect to the sample being perpendicular to the 
planes of the layers of the latter, 

means for modulating the amplitude of this excitation 
signal , 



25 
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a pyroelectric sensor in thermal contact with one side 
of the sample, this side being parallel to the planes 
of the said layers, 

means for measuring a photothermal voltage signal at 
5 the terminals of this pyroelectric sensor, 

characterised in that it comprises in addition: 

means for determining at least one numerical value 
representative of the amplitude of the said- 
photothermal signal and/or at least one numerical value 

10 representative of the phase difference between this 

photothermal signal and the excitation signal, 
calculation means in which is stored a two-dimensional 
axially symmetrical model of the thermal transfers that 
occur in the multilayer assembly consisting of the 

15 sample and the sensor, these calculation means enabling 

the value of the physical parameter corresponding to 
the measured .photothermal voltage signal to be 
determined . 

20 Other characteristics and advantages of the invention will 
also follow from the following description, which is purely 
illustrative and does not restrict the invention. This 
description should be read having regard to the 
accompanying drawings, in which: 

25 

Fig. 1 shows diagrammatical ly an arrangement for the 
implementation of the invention, 

Figs. 2 and 3 show diagrammat ically two possible 
3 0 analysis configurations, 



Fig. 4 illustrates the thermal model employed by the 
arrangement according to the invention, 



I I i 
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Fig. 5 illustrates an analysis configuration used to 
verify the correlation between the model proposed by 
the invention and the experimental results, 

5 

Figs. 6 and 7 are graphs on which are plotted, as a 
function of the parameter X shown in Fig. 5, on the one 
hand the amplitude arid the phase of the photothermal 
o .; voltage signal recorded with the arrangement of Fig. 5, 

10 and on the other hand the theoretical curves of the 

amplitude and phase of the voltage signal calculated 
from the model employed by the arrangement of the 
invention . 

15 The arrangement shown in Fig. 1 comprises an excitation 

source 1 of 1 mW power, emitting a beam F of monochromatic 
light onto a sample 2 to be analysed. At the output of 
this source 1 the beam F is modulated by means that may 
consist of a conventional mechanical cho£E.er, or even, as 

2 0 is shown in Fig. 1, of an electrooptical modulator 

controlled by a unit 3 that supplies a modulation signal 
thereto. Other luminous excitation sources are of course 
also possible. In particular, an excitation by a simple 
LED is sufficient. Other types of modulation are also 

2 5 possible, such as in particular a modulation by acoustical - 
optical techniques . 

A pyroelectric sensor 4 is mounted on one side of the 
sample 2. This pyroelectric sensor 4 is for example of the 
30 conventional thin film type comprising polyvinyl idene 

fluoride (PVDF) . The voltage signal at the terminals of 
this sensor 4 is passed to a preamplifier 5 whose output 
signal is transmitted, at the same time as a signal derived 
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from the modulation signal, to a Stanford SR-850 type 
synchronous detection amplifier 6. 



The voltages sampled by the amplifier 6 are passed to the 
5 calculation means 7, at the input of which they are 

digitised. The calculation means 7 process the amplitude 
and the phase of the received signals in order to determine 
the physical parameter of the sample that it is desired to 
measure . 

10 

Reference will now be made more specifically to Figs. 2 
and 3, which illustrate two possible analysis 
configurations . 

15 In the configuration of Fig. 2, the multilayer sample 2 
directly receives the modulated beam F # the pyroelectric 
sensor 4 being arranged on the side of the sample 2 
opposite to that which receives the luminous flux. 

20 In another possible configuration, termed "inverse 

configuration"; it is the pyroelectric sensor 4 that 
receives the beam F. Its side facing opposite the source 1 
carries a black layer 8 that absorbs the thermal energy of 
the beam F . 

2 5 

For a more detailed description of the pyroelectric sensors 
and arrangements with which they are associated, reference 
may advantageously be made to the following articles: 

"A versatile inverse photopyroelectric (IPPE) technique 

3 0 and instrument for real time observation of the 

condensation of water vapor in the atmosphere' 7 - Mihai 
CHIRTOC, Dane BICANIC, Valer TOSA - Rev. Sci. Instrum. 
62 (9), September 1991 pp. 2257-2261; 
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"Pyroelectric Sensors for the Photothermal Analysis of 
Condensed Phases" - H. COUFAL and A. MANDELIS, 
Ferroelectrics 1991, Vol. 118, pp. 379-409. 

5 Conventionally, the sensor 4 is composed of a flat 

pyroelectric film of surface area A and width L, metallised 
on both its sides (small thickness deposits) . This film is 
characterised by a spontaneous electrical polarisation P, 
which is a function of the temperature T. As a result of 

10 this polarisation phantom charges Q are formed on th^ 

surface, the surface density signal a of which is equal to 
the perpendicular component of the polarisation P N . These 
coupled charges create an internal electrical field which, 
due to the displacement of the free charges that it causes, 

15 leads to the neutralisation of the surface charge. The 

thermal order perturbations modify this surface equilibrium 
by acting on the polarisation. These thermal variations of 
the polarisation are characterised, at a given temperature, 
by the pyroelectric coefficient p (in cm^K" 1 ) . Any coupled 

20 charge variation in the dielectric as a result of a 

variation in the temperature may therefore be detected by 
measuring the potential difference that is established 
between the metal armatures. 

2 5 The voltage at the terminals of the pyroelectric sensor is 
given by: 




in which 



Q<r) »f T(r.t)d* 
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where s r is the characteristic relative permittivity of the 
dielectric (equal to about 12 for PVDF) . 

Reference will now be made to Fig. 4. In accordance with 
5 the invention, the sample 2 and the sensor 4 in contact 
with the latter are regarded as forming an axially 
symmetrical element with two-dimensional thermal transfers. 
The axis of symmetry is the axis of the excitation beam F, 
perpendicular to the planes of the layers. 

10 

The two-dimensional and axially symmetrical model of the 
thermal transfers that take place within the assembly of 
the sample/sensor system enables the temperature field in 
this multilayer structure to be determined and allows a 
15 calculation relationship to be established between the 

photopyroelectric voltage signal and the excitation signal 
passed to the sample. 

The sample that will now be described is assumed to be in 
20 good th ermal contact with a surface S of t he armature of 

the detector. The light absorption in this sample/detector 
system and the conversion, assumed to be non-radiative, of 
the light energy into heat increases the temperature of the 
pyroelectric film (in a region that is assumed to be • 
25 limited to the surface S) . If the excitation is modulated 
sinusoidally, the resultant temperature variation produces 
in turn a periodic variation of the charge, accumulated in 
the pyroelectric, displaced in phased by an angle (J) with 
respect to the excitation flux. 

30 

The two-dimensional and axially symmetrical model of the 
thermal transfers that take place within the assembly of 
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the .sample/detector system has in the first place enabled 
the temperature field in this multilayer structure to be 
determined, and then the photopyroelectric signal to be 



calculated . 



5 



The expression for the excitation flux is: 



If the excitation is a laser with a Gaussian- type profile, 
the illumination is of the form: 



where I 0 'is the power density at the centre of the laser 
beam, r 0 is the excitation radius for which the value of (j) 0 
is divided by e (I 0 = P / 7ir 0 2 ) - The beam is perpendicular 
to the surface; the reflected radiation flux I±" and 

15 transmitted radiation flux Ii + that are propagated parallel 
to the axis of the cylinder should be taken into account in 
each of the layers of the system. The system of two 
equations resulting from the transfer equation should then 
be resolved in the two specific directions. If i is the 

2 0 index of the layer and pi is the given optical absorption 
coefficient relative to this layer, then we have 
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the solutions of which are: 
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The diffusion of the light is ignored here. The conditions 
on the limits in z allow us to calculate the coefficients Pi 
and Mi, and consequently to deduce therefrom the expression 
for the sources of heat generated within each of the 
5 layers i : 

Since the radiation fluxes do not in this case have radial 
components, these expressions become in the case of a 
10 modulated pulsation source co : 

Ei(ra)-iPiMie'[^e"-€J«» 

It is assumed that the absorption of the energy is complete 
and that it is wholly converted into heat. It is also 
15 assumed that adiabatic conditions apply to the side walls. 
The solution of the equation for the heat in each of the 
layers now enables us to calculate fully the temperature 
field in the sample and in the pyroelectric . This equation 
relating to. a layer i takes the form: 




20 

where 

- Ki is the thermal conductivity of the layer i, 

- Ti is the temperature of the layer i, 

- ri is the density, 

25. - ci is the specific heat, 

- Gi and Ei are the terms of sources of heat generated in 
the layer i . 
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The temperature is then expressed as the sum of three 
components : 

TfrxA) - T„ b +T c (r.*)+T.tr.*>el- 

where 

5 T am b is the ambient temperature, 

T c (r,z) is the continuous component of the temperature, 
T a (r,z) is the alternative component at any point in the 
layer i . 

10 In order to determine the photopyroelectric signal , only 
the dynamic component of the temperature has to be taken 
into account. This has the following form: 

•»© 

where 

15 r c being the radius of the sample, 



bn is the nth zero of J±, 



20 



at is the thermal diffusivity of the medium in the layer i, 

and, finally, pi is the coefficient of optical absorption, of 
2 5 the layer i . 



The constants A i#n are determined by the conditions at the 
limits on z, namely: 
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mixed conditions and exchange coefficients at the front 
and rear side, 



conservation of the flux and equality of the 
temperatures at the interfaces. 



10 



The terms Ei >n and Gi >n are functions of, inter alia, 
coefficients I n obtained by resolving the terms of sources 
based on Bessel functions such as Jo : 



Mil, 



where 

&t -5- 

I.-it^ 

If the index corresponding to the layer of the pyroelectric 
15 is chosen for i, then one may deduce therefrom the 
expression for the pyroelectric voltage: 



V PPE (ci>)=. |q 



- 
















Pi. 



2j. 



11=1 ^n r c 



A(co) ei«* 



where 



X in ■ Aiafe - oJi - e - <*»l") + B 5n (e - c ^ 



20 q is the ratio of the sample/detector surfaces, 
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A (go) is the function transfer of the electrical circuit 
formed by the detector connected to a current -voltage 
preamplifier/converter . 

5 Since we are dealing with complex quantities in this case, 
this expression may accordingly be written in polar 
co-ordinates and the values of the amplitude and phase of 
the photopyroelectric signal may be derived therefrom. 

» 

10 The calculation means 7 store this model and are programmed 
so as to solve this equation in order to determine, from 
the amplitude and the phase that they receive, the physical 
parameter that it is desired to measure. 

15 In a possible variant the calculation means determine 

theoretical curves giving the phase or the amplitude as a 
function of the parameter that it is desired to measure, 
this parameter in turn being determined from the measured 
phase or amplitude by using these curves as a nomogram. 

20 

It should be noted that the proposed model enables the 
theoretical voltage at the terminals of the pyroelectric 
sensor to be determined from a knowledge of, on the one 
hand, the power and modulation of the exchange beam and, on 
25 the other hand, of the thermal diffusivity, absorption and 
thickness of each layer of the sample/sensor assembly. 

The arrangement according to the invention, is 
advantageously used t o follow c hanges in a physica l 
30 s tructure and in particular to follow sedimentation or 

I aye r ^g r o w t h phenomena. In this case the analysed sample 
is preferably, but not necessarily, modelled as a 
stratified bilayer. The wavelength of the excitation beam 
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is chosen so as to differentiate the two adjacent layers. 
The invention is used particularly advantageously to follow 
the sedimentation or coagulation of blood samples. 

5 In the case where an erythrocyte sedimentation is followed, 
the blood sample may to a first approximation be regarded 
as a stratified bilayer, consisting of a blood layer 
practically opaque to the excitation and a plasma layer 
semi-transparent to the wavelength of the excitation beam, 

10 the thickness of which increases with time. Knowing in 
particular the thermal conductivity of blood and plasma 
(respectively 0.58 and 0.56 W/m/K) , their density 
(respectively 1026 and 1108 kg/m 3 ) , their specific heat 
(4180 J/kg/K for both the plasma and blood) , it is possible 

15 to determine the thickness of the plasma layer at any point 
in time, as a function of the measurements of the amplitude 
and phase of the photothermal signal. The total thickness 
of the assembly formed by the sample to be analysed and the 
pyroelectric sensor is of the order of 1 mm. 

20 

In order to analyse a coagulation behaviour, the inverse 
configuration is preferably employed so that the sample is 
not optically disturbed. 

25 The invention may also advantageously be used for other 

applications. In particular it may be used to detect the 
appearance of bacteria or to follow the growth of cells, 
which phenomena cannot be handled by conventional 
photothermal techniques . 



Reference will now be made to Fig. 5, which shows the 
analysis configuration by means of which it has been 
verified that the model employed by the arrangement 
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according to the invention gives a good correlation with 
experimental values . 

The sample 2 is a water sample in which is immersed an 
5 aluminium film 9 with a blackened surface, which can be 
moved towards or away from the sensor 4 . The distance 
between the sensor 4 and the surface 8 is identified by X. 

The amplitude v and the phase (j) of the voltage signal 
10 measured at the sensor 4 were plotted for several values of 
the distance X on the graphs of Figs. 6 and 7. The 
curves Vc and <j)c of the amplitude and phase calculated by 
means of the model employed by the arrangement according to . 
the invention were also plotted on these graphs. These 
15 graphs show that the proposed model effectively enables one 
to obtain by calculation values that are in good 
correlation with the experimentally measured values. 

The invention has been described in the case in which the 
20 light beam incident on the sample is a sinusoidally 

modulated beam. However, the invention also applies to the 
case in which the sample is excited in a pseudo-random 
manner. Reference may advantageously be made in this 
connection to the French patent in the name of the 

2 5 applicant published under no. 2 6 63 74 5. 

It should be noted that one of the advantages of the 
invention is that it enables an analysis to be. carried out 
at a modul a t ion-~frie. que.ncjy: . b>el ow JI-Hz, which allows one to 

3 0 work with sample thicknesses of up to 1 mm. 
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Arrangement for the photothermal analysis of a 
material with a layered structure, for the 
determination of at least one physical parameter of 
one of these layers, such as its thickness, its 
thermal dif fusivity or its optical absorption, 
comprising : 

a source (1) for the emission of a luminous 
excitation signal in the direction of a sample (2) 
of the material to be analysed, the direction of 
incidence of this signal with respect to the 
sample being perpendicular to the planes of the 
layers of the latter, 

means (3) for modulating the amplitude of this 
excitation signal , 

a pyroelectric sensor (4) in thermal contact with 
one side of the sample, this side being parallel 
to the planes of the said layers, 

means (5, 6) for measuring a photothermal voltage 
signal at the terminals of this pyroelectric 
sensor, 

characterised in that it comprises in addition: 

means for determining at least one numerical value 
representative of the amplitude of the said 
photothermal signal and/or at least one numerical 
value representative of the phase difference 
between this photothermal signal and the 
excitation signal , 

calculation means (7) in which is stored a two- 
dimensional axially symmetrical model of the 
thermal transfers that occur in the multilayer 
assembly consisting of the sample (2) and the 
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sensor (4) , these calculation means enabling the 
value of the physical parameter corresponding to 
the measured photothermal voltage signal to be 
determined. 

Arrangement according to claim 1, characterised in 
that the means (5, 6) for measuring the photothermal 
voltage signal and for determining a numerical value 
representative of its amplitude and/or of its phase 
difference comprise means for storing the numerical 
values corresponding to photothermal , signals received 
at regular time intervals, in such a way that it is 
possible to follow evolving samples. 

Arrangement according to one or more of the preceding 
claims, characterised in that the calculation means 
model the flux of the excitation signal according to < 
Gaussian flux. 

Arrangement according to one of the preceding claims, 
characterised in that the model stored in the 
calculation means takes into account the transmitted 
and reflected radiation fluxes in each of the layers 
of the assembly formed by the sample and the sensor. 

Arrangement according to one of the preceding claims, 
for the analysis of a sedimentation or a layer growth 

Arrangement according to one or more of the preceding 
claims, characterised in that the calculation means 
model the sample as a stratified bilayer. 
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Arrangement according to one of claims 5 and 6, for 
the analysis of the sedimentation of a blood sample. 

Arrangement according to one of claims 5 and 6, for 
the analysis of the coagulation of a blood sample. 

Arrangement according to one of the preceding claims, 
characterised in that the excitation signal is 
modulated at a frequency of the order of or less than 
1 Hz. 
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